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EQUATIONS  OF  MOTION  OF  A  RIGID  PROJECTILE 

ABSTRACT 

The  basic  vector  equations  of  motion  of  a  rigid  body  are  described,  and 
an  aerodynamic  force^moment  system  is  developed.  The  vector  forces  and 
moments  due  to  gravity  and  rotation  of  the  earth  are  added  to  the  system. 
Forces  and  moments  produced  by  a  rocket  motor  are  included,  and  the  re¬ 
sulting  vector  differential  equations  of  motion  of  a  symmetric  rigid  projec¬ 
tile  are  derived. 

An  appendix  is  added,  which  describes  a  set  of  initial  conditions  imposed 
on  pertinent  parameters  for  rigid  body  trajectory  simulation. 
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Table  of  Symbols 


Term 

Definition 

Units 

A  (t) 

Axial  moment  of  inertia  at  time  t 

lb -ft2 

A  (t) 

Rate  of  change  of  A  at  time  t 

lb -ft2 
sec 

A 

e 

Area  of  jet  exit 

Az 

Azimuth  of  line  of  fire  (clockwise  from  north) 

deg 

B  (t) 

Transverse  moment  of  inertia  at  time  t 

lb -ft2 

2 

B  (t) 

Rate  of  change  of  B  at  time  t 

lb -ft 

sec 

CG 

Distance  from  nose  of  rocket  to  the  center  of 

mass 

ft 

d 

Reference  diameter  of  missile 

ft 

E 

Position  of  missile  with  respect  to  spherical 
earth  surface 

ft 

€ 

Fin  cant  angle 

rad 

f 

Deceleration  due  to  friction 

ft/  sec2 

g 

Acceleration  due  to  gravity 

ft/  sec2 

go 

Acceleration  due  to  gravity  (surface) 

ft/  sec2 

g 

Constant  used  in  the  conversion  of  thrust  (lbf) 

ft/  sec2 

c 

to  force  (lb~ft/sec  ) 

H 

Total  angular  momentum 

lb-ft2-rad/ 

h 

Angular  momentum  divided  by  B 

rad/ sec 

h 

Rate  of  change  of  h 

rad/ sec 

sec 
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Term 

\ 

\ 

MV 

XSP 

JST 

KA 

K 


K 


D 


K 


E 


K, 


K 


H 


K„ 


K 


M 


K, 


kt 

^F 


K 


XT 


e 

L 

M 


Definition 
h  at  end  of  launch  (tT  ) 

U  i-J 

h  at  end  of  launch  (t  ) 

3  Lj 

h  at  end  of  launch  (tT  ) 

-*  E 

Specific  impulse  per  unit  of  fuel  mass 

Total  impulse  at  standard  conditions 
Spin  damping  moment  coefficient 
Drag  coefficient 

Yaw  drag  coefficient 

Fin  cant  coefficient 

Magnus  force  coefficient 

Damping  moment  coefficient 

Lift  force  coefficient 

Overturning  coefficient 

Pitching  force  coefficient 

Magnus  moment  coefficient 

Magnus  cross  force  coefficient 

Magnus  cross  moment  coefficient 

Length  of  missile 

Latitude  of  launch  point 
Mach  number 


Units 
rad/ sec 

rad/ sec 

rad/ sec 

lbf-sec 

lbm 

lbf-sec 


l/  rad^ 


ft 

deg 
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Term 

Definition 

Units 

mf  (t) 

Mass  of  fuel  at  time  t 

lb 

m  (t) 

Mass  of  missile  at  time  t 

lb 

m  (t) 

Rate  of  change  of  mass  at  time  t 

lb/  sec 

N 

Axial  spin 

rad / sec 

P 

a 

Static  atmospheric  pressure 

lb/ft2 

P 

e 

Jet  pressure  at  nozzle  exit 

lb/ft2 

r 

Distance  between  center  of  earth  and  body 

ft 

ri 

Distance  from  CG  (t  )  to  CG  of  fuel 

33 

ft 

r2 

Radius  of  gyration  of  fuel  mass 

ft 

r 

e 

Distance  from  CG  of  missile  to  nozzle  exit 

ft 

rF 

Distance  from  tail  to  CG  of  fuel 

ft 

r 

s 

Radius  of  spin  rocket  action 

ft 

rt 

Distance  from  CG  of  missile  to  nozzle  throat 

ft 

r3(t) 

Distance  of  CG  shift  due  to  fuel  mass  change 
at  time  t 

ft 

t 

Time 

sec 

T* 

Effective  thrust 

lb 

Motor  burnout  time 

sec 

lBM 

Time  from  ignition  until  first  motion 

sec 

^MST 

Time  from  ignition  until  first  motion  at 
standard  condition 

sec 

^BST 

Motor  burnout  time  at  standard  condition 

sec 

tf 

Thrust  factor 

_ 
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m 

1  erm 

Definition 

Units 

fcL 

Time  at  end  of  launch 

sec 

tr  lt> 

Thrust  produced  by  motor  at  time  t 

Ibf 

T 

s 

Spin  rocket  thrust 

Ibf 

T 

Thrust 

Ibf 

u 

Velocity  of  missile  with  respect  to  ground 

ft/  sec 

u 

Acceleration  of  missile  with  respect  to  ground 

ft/  sec^ 

V 

Velocity  of  missile  with  respect  to  air 

it/ sec 

XCG  !t> 

Distance  from  nose  of  missile  to  center  of  mass 
at  time  t 

ft 

XCP 

Normal  force  center  of  fires  sure  from  nose 

ft 

X 

Position  of  missile  with  respect  to  ground 

ft 

X 

Unit  vector  along  longitudinal  axis  of  missile 

— 

X 

Rate  of  change  of 

rad/ se 

X  (t  ) 

JL 

x  at  end  of  launch  (t  ) 

L 

y 

Unit  vector  perpendicular  to  x 

-*■ 

y 

Rate  of  change  of  y 

rad/se 

5.V 

y  at  end  of  launch  (t  ) 

Ij 

2 

Unit  vector  perpendicular  to  both  x  and  y, 
z  =  (x  X  y) 

jngjgy  iai$t 

— 

z 

Rate  of  change  of  z 

z  (t  ) 

<*4*  Xj 

z  at  end  of  launch  (t,  ) 

rad/sei 
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Term 

Definition 

Units 

B 

a 

Yaw  of  missile 

deg 

7 

Angle  between  the  x  y  plane  and  the  angular 

deg 

r 

thrust  malalignment 

r 

Angle  between  the  x  y  plane  and  the  linear 

deg 

thrust  malalignment 

A 

Linear  thrust  malalignment 

ft 

6 

Angular  thrust  malalignment 

deg 

A 

**► 

Acceleration  due  to  rotation  of  earth 

ft/sec 

P 

Air  density  (varies  with  altitude) 

lb/ft3 

0 

Angle  of  elevation  of  launch 

deg 

Orientation  of  yaw 

deg 
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Introduction 


The  purpose  of  this  report  is  to  present  a  general  and  flexible  mathe¬ 
matical  model  for  rigid  body  trajectory  simulation. 

The  total  vector  force-moment  system  presented  is  based  on  a  collection 
of  the  ideas  of  several  authors.  In  a  basic  paper  published  in  1920  [Ref,  l]  , 
Fowler,,  Gallop,  Lock,  and  Richmond  developed  the  basic  aerodynamic  the~ 
ory  and  the  concise  vector  notation  of  the  present  report.  (See  also  Reference 
0,  and  unpublished  work  by  A.  S,  Galbraith).  Fowler’s  aerodynamic  hypoth¬ 
esis  contained  a  logical  contradiction,  as  was  pointed  out  by  Nielsen  and 
Synge  in  1946  [Ref.  2]  .  Kelley,  McShane,  and  Reno  in  their  text  on  exteri¬ 
or  ballistics  [Ref,  4]  extended  the  work  of  Nielsen  and  Synge  and  developed 
the  basic  aerodynamic  hypothesis  used  in  the  present  report. 

The  force-moment  equations  produced  by  a  rocket  motor  were  obtained 
by  collecting  the  basic  principles  expounded  in  two  texts  [Ref.  3,  6]  . 

The  resulting  equations  were  brought  into  logical  consistency  and  expressed 
in  the  vector  notation  of  the  present  report. 

The  usefulness  of  any  mathematical  model  is  determined  by  its  ability  to 
closely  match  the  results  of  physical  experiments  over  a  spectrum  of  test 
conditions*  The  model  presented  in  this  report  has  successfully  simulated 
free  flight  rocket  trajectories  over  a  wide  spectrum  of  test  conditions.  A 
limited  amount  of  simulation  has  been  performed  on  conventional  spin  stabi¬ 
lized  projectiles,  and  the  results  are  quite  encouraging. 
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BASIC  LAWS  OF  MOTION  OF  A  RIGID  BODY 


Consider  a  ground-fixed,  right-handed  coordinate  system,  whose  origin 
is  located  at  the  center  of  mass  of  a  rigid  body,  as  shown  in  Fig,  1  below. 


Assume  that  the  body  can  be  considered  a  solid  of  revolution,  and  assign 
to  the  axis  of  rotational  symmetry  a  unit  vector  x  in  the  chosen  coordinate 
system.  Since  the  rigid  body  is  to  represent  a  rocket  or  shell,  the  direction 
of  x^frorn  tail  to  nose  is  defined  as  positive. 

The  total  angular  momentum  of  the  body  can  now  be  expressed  as  the  sum 
of  two  vectors  in  the  ground-fixed  coordinate  system: 

(a)  The  angular  momentum  about 

(b)  The  total  angular  momentum  about  an  axis  perpendicular  to  x  . 
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The  angular  momentum  about  has  the  magnitude  AN,  where  A  is  the 
moment  of  inertia  of  the  body  about  Xj  and  N  is  the  axial  spin  or  angular 
velocity  about  xj  ^  ,  in  radians  per  second,  hence,  the  total  angular  momen¬ 
tum  about  x  can  be  represented  by  the  vector  AN  x^. 

The  total  angular  velocity  of  the  body  about  an  axis  perpendicular  to  x,is 
given  by  the  vector  (x^X  x),  where  the  superscript  dot  refers  to  differentia¬ 
tion  with  respect  to  time.  Since  the  body  possesses  rotational  symmetry, 
every  axis  through  the  center  of  mass  and  perpendicular  to  x  is  a  principal 
axis  of  inertia.  If  the  moment  of  inertia  of  the  body  about  any  transverse 
axis  is  B,  the  total  angular  momentum  about  an  axis  perpendicular  to  x  has 
the  vector  representation  B  (x^X  x). 

Let  H  denote  the  total  vector  angular  momentum  of  the  body.  The  vector 
representation  of  H  is: 

(1.  1)  H  =  ANx  +  B  (x  x  jy 

Let  h^=  H/B,  and  divide  both  sides  of  Equation  (1.  1)  by  B. 

AN 

(1.2)  h  =  ^x+  (xXx) 

Let  SM  be  the  sum  of  the  vector  applied  moments,  and  set  EM  equal  to 
the  vector  rate  of  change  of  angular  momentum. 

(1.  3)  h  =  ^  X  +  —  X+  (X  Xx}  =  s  (M/B) 


[  l]  In  this  report  a  positive  N  is  defined  as  rotation  which  would  cause  a 
right-hand  screw  to  advance  in  the  direction  of  x^. 
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In  addition  to  these  basic  equations,  two  additional  expressions  are  needed 
for  use  in  the  force^moment  system. 


(1.4)  (h  X  x)  =  x 

AN 

(1.5)  (h  .  x) 


Equation  (1.3)  is  the  basic  vector  differential  equation  of  angular  motion 
in  a  fixed  coordinate  system.  The  basic  equation  of  motion  of  the  center  of 
mass  is: 

(1.  6)  u  =  S  (F/m), 

where  2F  denotes  the  sum  of  the  vector  applied  forces,  m  is  the  mass  of 
the  body,  and  \j^is  the  vector  acceleration  of  the  center  of  mass  in  the  fixed 
coordinate  system. 

It  is  now  necessary  to  determine  the  forces  and  moments  which  comprise 
SF  and  2M,  In  the  next  section  the  basic  aerodynamic  forces  and  moments 
are  considered;  those  of  non-aerodynamic  origin  are  deferred  until  a  later 
section. 


AERODYNAMIC  FORCES  AND  MOMENTS 

A  force  or  moment  is  defined  to  be  aerodynamic  in  origin  if  it  is  produced 
by  interaction  of  a  rigid  body  and  its  atmospheric  medium.  The  aerodynamic 
forces  and  moments  presented  in  this  report  are  those  considered  necessary 
to  form  a  logically  consistent  mathematical  model  for  a  rigid  body  possessing 
rotational  symmetry. 
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DRAG 


Drag  is  historically  defined  as  resistance  to  forward  motion  of  a  projec¬ 
tile.  The  magnitude  of  drag  force  is  represented  in  classical  exterior  bal- 
2  2 

listics  as  pd  v  ,  where  p  is  the  density  of  the  atmospheric  medium, 
o 

d  is  the  reference  diameter  of  the  projectile,  v  is  the  forward  velocity  of  the 

center  of  mass  with  respect  to  the  air,  and  is  a  dimensionless  number 

o 

called  the  drag  coefficient.  If  the  forward  velocity  with  respect  to  air  is 
represented  by  the  vector  v,  then  the  direction  of  drag  is  —  v. 

Figure  2  illustrates  the  general  case  of  a  projectile  with  the  unit  vector 
pointing  along  the  axis  of  symmetry,  and  the  forward  velocity  represented 
by  v  .  The  angle  between  v  and  x,  denoted  by  a,  is  traditionally  called  the 
angle  of  yaw. 

Yaw  increases  drag  on  a  projectile  by  presenting  to  the  air  stream  an 

enlarged  cross-sectional  area.  The  effect  of  yaw  on  drag  is  accounted  for 

by  allowing  to  increase  with  yaw  squared.  The  yaw  drag  coefficient  is 

o 

denoted  by  ,  and  the  increase  in  the  drag  coefficient  due  to  yaw  is  then 

CL 

2  2 
given  by  a  ,  The  effective  drag  coefficient  is  (K^  +  a  ).  Based  on 

a  o  a 

the  above  definitions,  drag  force  is  represented  by  the  vector  equation: 

(2.  1)  DRAG  FORCE  =  -  pd2  (KD  +  KD  aZ)  v  v 

o  a 

SPIN  DAMPING  MOMENT 

The  spin  damping  moment  is  an  aerodynamic  moment  produced  by  vis¬ 
cous  friction  of  the  air  on  ths  surface  of  a  spinning  shell,  and  is  a  couple 
tending  to  destroy  axial  spin.  (See  Figure  3.  ) 
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The  magnitude  of  the  spin  damping  moment  is  given  by  pd  K  Nv,  where 

a 

is  a  dimensionless  number  called  the  spin  damping  moment  coefficient. 
Since  the  spin  damping  moment  opposes  axial  spin,  its  direction  is  given  by 
-  x.  The  vector  representation  of  the  spin  damping  moment  is  given  by: 


(2.  2)  SPIN  DAMPING  MOMENT  =  -  pd4  K  N  v  x 

v  A  •* 


Replacing  N  with  its  equivalent  from  Equation  (1.  5): 

4 

(2.  3)  SPIN  DAMPING  MOMENT  =  -  K  v  (h  •  x)  x 

A  A 


FIN  CANT  MOMENT 

If  a  projectile  is  symmetrically  finned,  and  the  fins  form  equal  non-zero 
angles  with  the  axis  of  symmetry,  a  pure  spin  generating  couple  will  be  pro¬ 
duced  during  flight.  The  fin  cant  angle,  e,  is  defined  as  positive  if  it  gener~ 
ates  a  positive  spin.  Figure  4  illustrates  the  fin  cant  moment  for  a  positive 
fin  cant  angle. 

3  2 

The  magnitude  of  the  fin  cant  moment  is  given  by  pd  K  e  v  ,  where  K 

E  E 

is  the  dimensionless  fin  cant  moment  coefficient,  and  e  is  the  fin  cant  angle. 
The  direction  of  the  moment  is  x  for  positive  values  of  e.  Accordingly,  the 
vector  representation  of  the  fin  cant  moment  is: 

(2.  4)  FIN  CANT  MOMENT  =  pd3  K  e  v2  x 
LIFT  FORCE 

Aerodynamic  lift  is  created  by  an  asymmetric  air  flow  over  a  yawed  body. 

2  2 

The  magnitude  of  the  lift  force  is  given  by  pd  K  v  sin  a,  where  K  is  the 

E  E 

dimensionless  lift  force  coefficient.  The  lift  force  acts  in  the  plane  of  yaw 
and  is  perpendicular  to  the  direction  of  motion  of  the  projectile. 
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Figure  5  illustrates  the  vector  representation  of  the  lift  force*  Consider 

2 

the  vector  [  v^X  (x^X  v)]  .  This  vector  has  the  magnitude  v  sin  a,  and  is 
perpendicular  to  v^in  the  plane  containing  v  and  x^  The  vector  representa^ 
tion  of  lift  force  is  then  given  by: 

(2.  5)  LIFT  FORCE  =  pd2  [  v  X  (x  X  v)] 

Expanding  the  triple  vector  product,  and  rewriting: 

(2.  6)  LIFT  FORCE  =  pd2  K  [  v2  x  -  (v  .  x)  v] 

The  vector  sum  of  the  lift  and  drag  forces  can  be  expressed  as  a  total  re¬ 
sistance  vector*  If  the  components  of  vector  resistance  are  resolved  paral¬ 
lel  and  perpendicular  to  x  instead  of  v,  the  axial  drag  force  and  normal  force 
result.  The  normal  force  is  the  component  of  the  vector  resistance  that  pro¬ 
duces  the  overturning  moment. 

OVERTURNING  MOMENT 

If  the  line  of  action  of  the  aerodynamic  normal  force  does  not  pass  through 

the  center  of  mass  of  the  projectile,  an  overturning  moment  will  be  produced. 

The  magnitude  of  the  overturning  moment  due  to  normal  force  is  given  by 
3  Z 

pd  K  ,  v  sin  a,  where  K  ,  is  the  dimensionless  overturning  moment  coef- 
M  M 

ficient.  The  overturning  moment  is  perpendicular  to  the  plane  of  yaw,  or  to 
both  v  and  x^  as  is  illustrated  in  Figure  6.  The  vector  (v^X  x)  has  the  proper 
direction,  and  in  magnitude  is  equal  to  v  sin  c.  The  overturning  moment 
vector  is  then  given  by: 

(2.  7)  OVERTURNING  MOMENT  =  pd3  v  (v  X  x) 
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If  the  line  of  action  of  normal  force  intersects  the  axis  of  symmetry  at  a 
point  behind  the  center  of  mass,  the  overturning  moment  acts  as  a  restoring 
moment.  This  situation  is  met  by  allowing  K  to  be  negative.  The  addition 
of  fins  at  the  rear  of  a  body  of  revolution  moves  the  center  of  pressure  of  the 
normal  force  to  the  rear,  and  usually  results  in  a  negative  K.^,  or  restoring 
moment. 

MAGNUS  FORCE 

The  Magnus  force  arises  from  the  interaction  of  the  air  stream  and  the 

boundary  layer  of  a  yawed  spinning  body.  The  magnitude  of  the  Magnus  force 
3 

is  given  by  pd  K  N  v  sin  a,  where  K  is  the  dimensionless  Magnus  force 

X  X 

coefficient.  The  direction  of  the  Magnus  force  is  perpendicular  to  the  plane 
of  yaw  and  is  represented  by  the  vector  (x^X  v)  for  positive  values  of  N,  as 
illustrated  in  Figure  7.  The  Magnus  force  is  represented  by  the  vector: 

(2.  8)  MAGNUS  FORCE  =  pd3  Kp  N  (x  X  v) 

Replacing  N  with  its  equivalent  from  Equation  (1.  5): 

(2.  9)  MAGNUS  FORCE  =  (h  .  x)  (x  X  v) 

MAGNUS  MOMENT 

If  the  line  of  action  of  the  Magnus  force  does  not  pass  through  the  center 

of  mass  of  the  projectile,  a  Magnus  moment  will  be  produced.  The  magni- 

4 

tude  of  the  Magnus  moment  is  given  by  pd  K^,  N  v  sin  a,  where  K^,  is  the 
dimensionless  Magnus  moment  coefficient.  The  Magnus  moment  lies  in  the 
plane  of  yaw,  and  is  perpendicular  to  x.  The  vector  [  3^X  (x^X  V)]  has  the 
magnitude  v  sin  a  and  has  the  proper  direction,  as  illustrated  in  Figure  8. 

The  Magnus  moment  is  represented  by  the  vector: 

(2.  10)  MAGNUS  MOMENT  =  pd4  KT  N  [  x  X  (x  X  v)] 
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Again  replacing  N  with  its  vector  equivalent,  and  expanding  the  vector 
triple  product: 

4 

(2.  11)  MAGNUS  MOMENT  =  K  (h  *  x)[(v  *  x )  x  -  v  ] 

PITCHING  FORCE 

The  pitching  force  is  a  force  opposing  any  change  in  the  direction  of  the 
longitudinal  axis  of  a  projectile*  If  the  unit  vector  x  is  in  motion,  its  rate 
of  change,  is  in  the  plane  of  motion,  and  is  perpendicular  to  x^.  If  |^x| 
represents  the  scalar  magnitude  of  x,  the  magnitude  of  the  pitching  force  is 
pd  K  v|x|,  where  K  is  the  dimensionless  pitching  force  coefficient*  The 
direction  of  the  pitching  force  is  parallel  to  x  and  oppositely  sensed,  as  is 
illustrated  in  Figure  9*  The  vector  representation  of  the  pitching  force 
follows : 

3 

(2.  12)  PITCHING  FORCE  =  -  pd  K  v  x 

o  ^ 

Replacing  x  with  its  equivalent  from  Equation  (1.  4): 

(2.  13)  PITCHING  FORCE  -  -  pd3  K„  v  (h  X  x) 

DAMPING  MOMENT 

The  damping  moment  is  a  moment  opposing  angular  velocity  of  the  long¬ 
itudinal  axis  of  a  projectile,  and  is  the  aerodynamic  moment  associated  with 
pitching  force.  If  the  unit  vector  x_  is  in  motion,  its  vector  rate  of  change  is 
x,  and  the  transverse  angular  velocity  vector  is  X  x).  The  magnitude  of 
the  damping  moment  is  pd4  K  v  |x|,  where  K  is  the  dimensionless  damp- 
ing  moment  coefficient*  The  direction  of  the  damping  moment  is  parallel  to 
(x  X  x)  and  oppositely  sensed,  as  is  illustrated  in  Figure  10*  The  damping 
moment  is  represented  by  the  vector: 

(2.  14)  DAMPING  MOMENT  =  -  pd4  K „v  (x  X  x) 

ri  ^ 
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Replacing  x  with  its  vector  equivalent,  and  expanding  the  resulting  triple 
vector  product: 

(2.15)  DAMPING  MOMENT  =  -  pel4  KR  v  [h  -  (h  .  x)  x] 

MAGNUS  CROSS  FORCE 

The  Magnus  cross  force  is  a  magnus  force  arising  from  a  transverse 
angular  velocity  of  the  shell  axis*  The  magnitude  of  the  Magnus  cross  force 
is  pd^  N  |  x|  ,  where  K^..^  is  the  dimensionless  Magnus  cross  force  co¬ 

efficient.  The  direction  of  the  Magnus  cross  force  is  parallel  to  the  vector 
(x^X  x),  and  has  the  same  sense  for  positive  values  of  N,  as  shown  in  Figure 
11.  The  Magnus  cross  force  is  represented  by  the  vector: 

(2.16)  MAGNUS  CROSS  FORCE  =  pd4  Kv„  N  (x  X  x) 

yv  jc  ~L* 

Replacing  N  and  x^with  their  respective  vector  equivalents,  and  expand^ 
ing  the  resulting  triple  vector  product: 

4 

(2.  17)  MAGNUS  CROSS  FORCE  =  Kv„  (h  .  x)  [  h  -  (h  .  x)  xl 

MAGNUS  CROSS  MOMENT 

The  Magnus  cross  moment  is  the  aerodynamic  moment  associated  with 

the  Magnus  cross  force.  The  magnitude  of  the  Magnus  cross  moment  is 
5  ,  # , 

pd  K  N  |x|,  where  K  is  the  dimensionless  Magnus  cross  moment  co- 
X.  1  *  X 1 

efficient.  Since  the  direction  of  the  Magnus  cross  moment  is  parallel  to  x^ 
and  oppositely  sensed,  as  is  shown  in  Figure  12,  the  Magnus  cross  moment 
is  represented  by  the  vector: 

(2.  18)  MAGNUS  CROSS  MOMENT  =  -  pd§  KXT  N  x 
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Replacing  N  and  x  with  their  respective  vector  equivalents: 

^  5 

(2.  19)  MAGNUS  CROSS  MOMENT  =  -  pd  -  B  K  (h  .  x)  (h  x  x) 

The  dimensionless  aerodynamic  coefficients  are  functions  of  many  dimen¬ 
sionless  power  products,  including  the  dimensionless  shape  parameters, 
Reynolds  number  and  Mach  number.  Aerodynamic  coefficients  are  defined 
with  reference  to  a  specific  set  of  shape  parameters,  and  may  be  expressed 
as  functions  of  Mach  number.  Additional  aerodynamic  coefficients  may  be 
included  to  account  for  the  variation  of  aerodynamic  forces  and  moments  with 
yaw. 
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MAGNUS  MOMENT 


GRAVITY  AND  ROTATIONAL  ACCELERATIONS 


The  acceleration  due  to  gravity  is  caused  by  the  force  of  gravitational 
attraction  between  a  body  and  the  earth*  The  magnitude  of  gravitational  force 
is  proportional  to  the  mass  of  the  body,  and  inversely  proportional  to  the 
squared  distance  between  the  centers  of  mass  of  the  body  and  of  the  earth. 

The  line  of  action  is  directed  from  the  mass  center  of  the  body  to  the  center 
of  mass  of  the  earth. 


It  is  assumed  that  a  spherical  earth  is  a  sufficiently  good  approximation 
for  the  purpose  of  trajectory  simulation.  The  position  of  the  missile  center 
of  mass  is  specified  by  the  vector  X  in  the  ground  fixed  coordinate  system* 
The  force  of  gravity  is  approximated  by  the  vector: 


(3-D  g  =  “  g 


R~ 

o  3 
r 


(X2  +  R) 


X. 


r  2  , 2  Z,  1/2 

[X^  +  (X  +  R)  +  X*  ]  ' 


where  g  -  value  of  gravitional  acceleration  at  point  of  launch 
R  =  radius  of  the  earth  at  point  of  launch 
r  =  distance  between  center  of  earth  and  body 

X  ,  X  ,  X  =  components  of  X  in  the  ground  fixed  coordinate  system. 

L  L*  3 
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In  the  derivation  of  the  force-moment  system,  all  vectors  have  been  ref¬ 
erenced  to  the  ground  fixed  coordinate  system.  Since  the  earth  is  spinning 
about  an  axis  passing  through  its  center  of  mass,  the  acceleration  produced 
by  rotation  of  the  earth  must  be  added  to  the  force  equation.  The  accelera¬ 
tion  due  to  rotation  of  the  earth  is  represented  by  the  vector: 

(3.2)  A  = 

The  \'s  are  defined  by  the  following  equations: 

\  =  2  8  cos  L  sin  A 

1  z 

(3.  3)  \  =28  sin  L 

M 

=  2  ft  cos  L  cos  A 
3  z 

where  S2  =  Angular  velocity  of  the  earth  (radians/ sec) 

L  =  Latitude  of  launch  point 

A  =  Azimuth  of  fire  measured  clockwise  from  North 
z 

u^ ,  u^,  u^  =  Components  of  the  vector  velocity  of  the  missile  with 
respect  to  ground 


' 

1 

— 

u 

1  Z 

2  3 

u, 

+ 

\  u 

1  1 

3  3 

\  u 

2  1 

3  2 

*  The  \’s  given  are  for  the  northern  hemisphere.  For  the  southern  hemi¬ 
sphere  replace  L  by  (-L)* 
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JET  FORCES  AND  MOMENTS 


The  simulation  of  free  flight  rocket  trajectories  necessitates  the  inclusion 
in  the  simulator  of  the  forces  and  moments  produced  by  the  rocket  motor. 

The  rocket  is  assumed  to  be  a  rigid  body,  and  all  jet  forces  are  assumed  to 
originate  at  the  nozzle. 


THRUST 

The  thrust  of  the  rocket  motor  is  a  vector  assumed  rigidly  attached  to 
the  rocket,  and  is  allowed  to  vary  in  magnitude  as  a  function  of  time. 

Consider  two  thrust“time  functions  as  illustrated  in  Fig.  13  below. 
THRUST 


Figure  13, 
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The  relevant  symbols  are  defined  below: 

A  =  Area  of  jet  exit 
e 

F[  ]  indicates  na  function  ofM 
I  -  Specific  impulse  per  unit  of  fuel  mass 

S  -L 

I  =  Total  impulse  at  standard  conditions 

m.  =  Mass  of  fuel 
f 

P  =  Static  atmospheric  pressure 
a 

P  -  Jet  pressure  at  nozzle  exit 
e 

T*  =  Effective  thrust 

t  =  Motor  burnout  time 
B 

*BM  =  ign^ion 

^BMST  =  ^me  °i  ignition  at  standard  conditions 

t  =  Motor  burnout  time  at  standard  conditions 
B  ST 

Tp,  =  Thrust  factor 

T  (t)  ~  Thrust  produced  by  motor  at  time  t 

XV 

T  =  Thrust 


Let: 

(4.  1)  T  =  FI 


.  f  *B ST  ~  ^BMST 

”  BM  l  t  -  t 

\  B  BM 


+  Smst^J 


The  thrust  of  the  rocket  motor  as  a  function  of  time  may  be  linearly 
transformed,  maintaining  the  same  total  impulse,  by  the  following  equation: 


(4.2)  Tr  =  TFmf 


(rSP  \  /  *~BST  ~  ^MST  \ 

*ST  j  \  *B  ~  *BM  / 
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The  force  due  to  jet  pressure  at  the  rocket  nozzle  is  given  by  P  A  .  The 

e  e 

static  atmospheric  pressure  produces  a  force  equal  to  -P  A  .  The  total 

a  e 

thrust  contributed  by  pressure  is  then  given  by  (P  —  P^)  A^. 

The  effective  thrust  of  the  rocket  motor  is  the  algebraic  sum  of  the  com¬ 
ponents  considered,  and  is  given  by: 


I 


(4.  3)  T*  =  T  m 


SP 


F  f  \  I 


ST 


*BST  _  *BMST 

t  -  t 
B  BM 


+  (P  -  P  )  A 
e  a  e 


The  velocity-time  history  of  a  rocket  during  the  burning  period  is  repro¬ 
duced  in  the  simulator  by  linearly  transforming  the  standard  thrust  curve  in 

time,  and  by  iteratively  determining  the  value  of  T  required  to  match  the 

r 

observed  burnout  velocity- 


s 


I 
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Figure  14  is  a  diagram  of  the  thrust  vector  in  the  missile  fixed  coordi¬ 


nate  system* 


y 


Figure  14* 


The  point  0^  represents  the  position  of  the  center  of  mass,  and  0^  repre¬ 
sents  the  position  of  the  throat  of  the  rocket  nozzle*  The  distance  from  0^  to 
0  is  r  ,  and  is  defined  as  a  positive  number. 

Ct  x 

Two  types  of  rocket  thrust  malalignments  are  considered  in  this  report. 
The  first  is  historically  defined  as  linear  malalignment,  in  which  the  nozzle 
is  displaced  from  the  axis  of  symmetry  of  the  rocket.  The  amount  of  dis¬ 
placement  is  denoted  by  A,  and  the  direction  by  an  angle  F,  measured  clock¬ 
wise  from  the  plane  containing  x  and  y. 
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The  second  type  is  defined  as  angular  malalignments  in  which  the  direc¬ 
tion  of  effective  thrust  forms  a  fixed  angle  with  the  axis  of  symmetry.  This 
direction  in  the  missile  fixed  coordinate  system  is  specified  by  two  angles; 
the  magnitude  of  angular  malalignments  denoted  by  6,  and  the  direction  angle 
y  *  measured  clockwise  from  the  plane  containing  x  and  y. 

The  vector  thrust  force  can  be  resolved  into  components  along  the  missile 
fixed  coordinate  system,  and  the  following  expression  results: 

(4.  4)  THRUST  FORCE  as  g^  T*  (cos  6  x^+  cos  y  sin  6  y^+  sin  y  sin  <5  z) 

Consider  a  vector  normal  to  the  thrust  vector,  passing  through  0^,  and 
directed  from  0^  toward  the  vector  thrust.  For  small  values  of  the  angle  6, 
this  vector  is  approximated  by: 

[  (A  cos  E  +  rj  cos  7  sin  6)  y  +  (A  sin  V  +  r  sin  y  sin  6)  z] 

The  vector  thrust  moment  is  the  cross  product  of  this  vector  onto  the 
vector  thrust  force.  If  the  indicated  expansion  is  performed,  the  following 
equation  results: 

(4.  5)  THRUST  MOMENT  =  g  T*  {[(A  cos  T  +  r  cos  7  sin  6)  sin  7 

c  t 

-  (A  sin  T  4  r  sin  y  sin  6)  cos  y  ]  sin  6  x  +  (A  sin  F  4  r  sin  y  sin  6)  cos  6  y 

t  t 

-  (A  cos  r  +  r  cos  7  sin  6)  cos  6  zj 

Equations  (4.  3),  (4.  4)  and  (4.  5)  represent  the  total  vector  force  and  mo¬ 
ment  due  to  the  thrust  of  the  rocket  motor. 
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JET  DAMPING 


Consider  a  burning  rocket  which  has  acquired  a  transverse  angular  veloc” 

ity.  In  addition  to  a  loss  in  mass,  the  rocket  is  also  losing  transverse 

angular  momentum*  It  is  assumed  that  the  jet  stream  passes  from  subsonic 

to  supersonic  velocity  at  a  point  midway  between  the  throat  and  the  exit  of 

the  nozzle*  The  nozzle  is  then  removing  transverse  angular  momentum 

from  the  system  at  the  rate  m  r  r  IxL  where  r  is  the  distance  from  the 

e  t  e 

rocket  center  of  mass  to  the  nozzle  exit.  The  effect  is  analogous  to  a  damp” 
ing  force  and  moment,  and  the  terms  are  therefore  called  jet  damping  force 
and  moment. 

B 

The  jet  damping  force  has  the  magnitude  (^ ”  -  m  r  )  |xj,  and  the  direc- 

rt  6  ^ 

tion  is  parallel  to  3^.  The  jet  damping  force  is  represented  in  Figure  15. 

(4.  6)  JET  DAMPING  FORCE  =  (~ - m  r  )  x 

r  e  ^ 

t 

Replacing  x  with  its  equivalent  from  Equation  (1.  4): 

(4.  7)  JET  DAMPING  FORCE  =  (—  -  m  r  )  (h  X  x) 

r  e 

t 

Figure  16  illustrates  the  jet  damping  moment  associated  with  the  jet  damp¬ 
ing  force.  The  magnitude  of  the  jet  damping  force  is  given  by  (B  -  m  r^  r^) 
j^c| .  It  is  parallel  to  the  vector  X  and  oppositely  sensed. 

(4.  8)  JET  DAMPING  MOMENT  *  -  (B  -  m  r  r  )  (x  X  x ) 

e  t •*  ■* 

Replacing  x^  with  its  vector  equivalent,  and  expanding  the  resulting  triple 
vector  product: 

(4.  9)  JET  DAMPING  MOMENT  =  -  (B  -  m  r  r  )  [h  -  (h  *  x)  x  ] 

e  t  ^ 
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SPIN  ROCKET  MOMENT 


Spin  rockets  are  utilized  for  imparting  spin  to  a  ballistic  rocket  after  it 
has  left  the  launcher*  If  a  symmetric  bank  of  spin  rockets  is  mounted  per¬ 
pendicular  to  the  longitudinal  axis  of  the  rocket,  the  result  is  a  pure  spin 

generating  couple.  The  total  thrust  of  the  spin  rockets  is  denoted  by  T  ,  and 

s 

the  radius  of  spin  rocket  action  by  r  ,  The  magnitude  of  the  spin  rocket  mo¬ 
ment  is  given  by  g  r  T  ,  where  g  is  a  constant  used  in  the  conversion  of 

c  s  s  2  c 

thrust  (lbf)  to  force  (lb-ft/sec  ).  The  direction  of  the  spin  rocket  moment  is 
x^  for  a  positive  spin  (Fig*  17).  The  vector  representation  of  the  spin  rocket 
moment  is: 

(4.  10)  SPIN  ROCKET  MOMENT  =  g  r  T  x 

c  s  s 

The  thrust  of  the  spin  rockets  is  linearly  transformed  in  time  by  the  same 
method  previously  described  for  the  rocket  motor* 
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EQUATIONS  OF  MOTION 

In  the  previous  four  sections,  the  forces  and  moments  acting  on  a  shell 
or  rocket  have  been  derived.  If  these  forces  and  moments  are  summed  and 
substituted  in  equations  (1.6)  and  (1.  3),  the  following  vector  differential 
equations  of  motion  result: 


(5.1)  u  =  “  (K  +  K  a2)  v  v  +  K  [  v^x  -  (v  .  x)  v  ] 

o  a 


2 

m  L 


3  3 

-  K  v  (h  x  x)  +  JL  K  (h  .  x)  (x  X  v) 

m  s  m  A  F  -* 


4 

+  £^~  a"  kxf  (V  $  til"  Ql*  2)  sJ  +  % 


m 


gcT* 

+  A  + -  (cos  6  x  +  cos  7  sin  6  y  +  sin  7  sin  6  z) 

-+  m  -+  — - 


+  (^~  —  m  r  )  (h  X  x) 

m  r  e 
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(5.2)  4  = 


j 3  ,4 

•  K  v  (v  X  x)  ~  K  v  [h  -  (h  .  x)  xl 

B  M  -*  -*  B  H  l~+  -*•  -*  -r 


d4  d3  2 

^  K  v  (h  .  x)  x  +  —  K  e  v  x 

A  A  b  E  •+ 


+  £^_  K  (h  .  x)  [(v  .  xj  x  -  vj 


,5  g  r  T 

—  —  K  (h  -  x)  (h  X  x)  +  - — 

A  XT  V  B 


gc  T* 

+  — {[  (A  cos  F  +  cos  7  sin  6)  sin  y 

-  (A  sin  r  +  rt  sin  y  sin  6)  cos  y  ]  sin  6  x^ 


+  (A  sin.  F  +  rf  sin  7  sin  6)  cos  6  y  (A  cos  T 
+  cos  7  sin  6)  cos  6  zj 


B  _  rrTr  r 

~  ( - - — - — ~)  [h  -  (h  .  x)  x] 


In  addition  to  these  two  basic  differential  equations,  the  equations  relat~ 
ing  the  (x,  y,  z)  coordinate  system  and  the  ground  fixed  coordinate  system 
are  needed.  Equation.  (1,4)  gives  (h  x3)-  The  rates  of  change  of  y^and 
z  in  the  ground  fixed  system  are  given  by: 
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B  —  A 

(5.3)  y  =  (hXy)  +  (— )  (h  .  x)  z 

(5.4)  ^  =  (hx  z)  »  (^“)  (h.  x)  y 

The  position  of  the  missile  center  of  mass  in  the  ground  fixed  coordinate 
system  is  represented  by  the  vector  X  ,  and  is  given  by  the  following  integral 
equation: 

(5.5)  X  =  f  u  dt 

^  o 

The  missile  center  of  mass  position  with  respect  to  the  spherical  earth 
surface  is  denoted  by  E.  The  vector  E  is  related  to  the  radius  of  the  earth, 
R,  and  to  X  by  the  following: 
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The  physical  properties  consisting  of  mass?  center  of  mass,  arid  moments 
of  inertia  must  be  considered  functions  of  time  for  rocket  trajectory  Simula-  i 

tion.  The  rate  of  change  of  mass  with  respect  to  time  is  given  at  any  instant 
by  the  ratio  of  motor  thrust  to  motor  specific  impulse.  The  change  in  center  *  * 

of  mass  and  axial  moment  of  inertia  are  represented  by  linear  functions  of 
time,  and  the  transverse  moment  of  inertia  is  found  by  applying  the  parallel 
axis  theorem  to  the  changing  mass.  The  following  equations  and  definitions 
summarize  the  method  used  to  account  for  time  variation  of  the  rocket’s  phys¬ 
ical  properties, 

xR(t) 

(5.7)  m(t)  =  -r— 

SP 
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(5.  13) 


mf  (t)  r  i  =  m  (t)  r3  (t) 

Where  r  =  Distance  from  CG  (t  )  to  CG  of  fuel 

I 

r2  ~  -^ac^us  gyration  of  fuel  mass 

r^(t)  =  Distance  of  CG  shift  due  to  fuel  mass  change  at  time  t 

X  (t)  =  Distance  from  nose  of  missile  to  center  of  mass  at 
GCj 

time  t 

YAW  AND  ORIENTATION  OF  YAW 

The  yaw  of  a  projectile,  denoted  by  a,  is  the  included  angle  between  the 
vectors  ^and  v^,  and  is  always  defined  as  a  positive  quantity*  Since  the  dot 
product  of  two  vectors  equals  the  product  of  their  magnitudes  multiplied  by 
the  cosine  of  their  included  angle,  the  magnitude  of  yaw  is  given  by: 

(5*  14)  a  =  Cos 


The  orientation  of  yaw  is  the  angle  between  the  plane  containing  both  v^ 
and  x  and  a  vertical  plane  containing  v^,  and  is  measured  clockwise  from  the 
vertical  plane.  Yaw  orientation  is  denoted  by  ^  >  and  is  given  by  the  following 
expre  s  sion: 

.*  2,1/2  „  v/2~ 

*3  (v  ~  v3  >  -  v3  ~  *3  > 

v  sin  a 


(5.  15)  i|>  =  Sin 
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4 j  becomes  indeterminate  when  a  is  zero  in  Equation  (5.  15).  If  4*  is  de¬ 
fined  as  zero  when  a  equals  zero,  Equation  (5,  15)  becomes  continuous  for 
all  values  of  a. 

F. 

Robert  F.  Lieske 

Robert  L.  McCoy 
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APPENDIX 


INITIAL  CONDITIONS. 

The  launching  phase  of  a  burning  rocket  is  characterized  by  a  period  of 
constrained  motion  prior  to  launching  into  free  flight.  It  is  assumed  that 
both  the  rocket  and  the  launcher  are  rigid  structures.  It  is  further  assumed 
that  the  rocket  can  acquire  no  angular  momentum  while  on  the  launcher,  and 
that  the  center  of  mass  is  constrained  to  move  in  a  straight  line  on  the  launch™ 
ing  rail. 


The  forces  permitted  to  act  on  the  rocket  during  launching  are  thrust, 
drag,  gravity,  and  rocket-launcher  friction.  If  0  is  the  angle  of  elevation  of 
the  launcher,  and  F  is  the  deceleration  due  to  rocket-launcher  friction,  the 
following  equations  describe  the  motion  of  the  rocket  center  of  mass  during 
the  launching  phase: 


(A.  1)  u 


gcT* 

m 

pd2 

COS  0  “  - 
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kdv2 

cos  0 

-  g  sin 
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^  2 

0  COS  0  -  F  COS  0 
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r  m 

1  "V 
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sin  0 

-  g  sin 
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0  -  F  cos  0  sin  0 

0 
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If  the  projectile  is  launched  from  a  gun  tube,  the  launching  phase  may  be 
omitted,  and  an  initial  velocity  is  assumed  at  end  of  launch. 


At  time  tT  ,  the  initial  conditions  of  axial  spin  and  transverse  angular 
momentum  are  assigned.  The  direction  cosines  of  the  missile  axis  of  sym¬ 
metry  are  given  by  the  components  of  3^.  The  unit  vector  y  is  defined  per¬ 
pendicular  to  3^,  and  lies  initially  in  the  vertical  plane  containing  3^.  The 
unit  vector  ^is  then  given  by  (x^X  y).  The  initial  values  of  3^,  y ,  and  ^are 
given  by  the  following  equations: 


(A.  2) 


(A.  3) 


(A.  4) 


COS  0 

sin  0 
0 


h?  +  h«  are  the  end  of  launch  components  of  l^with  respect  to  the  shell 
1  'l 


fixed  axis  system.  The  initial  h  used  for  simulation  must  include  the  effects 

of  initial  angle  of  elevation  of  the  trajectory  and  initial  axial  spin.  The  com- 

AN 

ponent  of  h^due  to  axial  spin  is  x.  The  initial  value  of  l^at  end  of  launch 

ig  given  by: 


(A.  5) 


V* 
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Values  of  h 


<V 


are  chosen  so  as  to  allow  the  simulator  to  match  the  ob¬ 


served  initial  angular  momentum  of  the  projectile. 
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